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The application of the Lorentz integral transform (LIT) method to photon scattering off nuclei is 
presented in general. As an example, elastic photon scattering off the deuteron in the unretarded 
dipole approximation is considered using the LIT method. The inversion of the integral transform is 
discussed in detail paying particular attention to the high-energy contributions in the resonance term. 
The obtained El-polarizabilities are compared to results from the literature. The corresponding 
theoretical cross section is confronted with experimental results confirming, as already known from 
previous studies, that the El-contribution is the most important one at lower energies. 

PACS numbers: 21.45.-v, 21. 45. Be, 25.20.De 

I. INTRODUCTION 



^SJ ' The study of electromagnetic reactions in photoabsorption and photon scattering on nuclei is an excellent tool 

to investigate nuclear structure. In addition, it can also lead to valuable insights into the properties of the nuclear 
constituents, the nucleons, like for example, electric and magnetic polarizabilities. In this context photon scattering 
experiments are a particularly interesting source of information on off-shell properties. On the other hand genuine mi- 
croscopic calculations of photon scattering cross sections are rather complicated since the complete nuclear excitation 
spectrum has to be taken into account. Thus it hardly comes as a surprise that in the past such theoretical efforts 
were mainly concentrated on the two-nucleon system. A first realistic calculation of deuteron photon scattering has 
been carried out in [1|. In the last decade quite a few theoretical investigations have been performed [2-Q among 
them also calculations based on chiral effective field theory. 

Considering a more complex A-nucleon system, as mentioned, one needs to have under control the corresponding 
A-body continuum. Today this could in principle be realized for the three-nucleon system, but many-body calculations 
of photon scattering cross sections for systems with A > 3 are presently out of reach. On the other hand, there is 
^vj I a particular interest in ^Li photon scattering recently, and data have already been taken at HIGS [3]. Fortunately, 
fT^ . the problem of calculating the correct many-body continuum wave function can be avoided by application of the 
PsJ ' Lorentz integral transform (LIT) method [^, Q. In fact, the LIT approach reduces a scattering state problem to 
a simpler bound-state like problem, which then can be solved with techniques that usually are applied for bound 
(^ • states. This leads to an enormous reduction of the complexity of the calculation, e.g., LIT calculations of the total 
photoabsorption cross sections have even been made for the six- and seven-body nuclei taking into account all possible 
break-up channels and with full consideration of the final state interaction [iO,lll|- In view of the fact, that till now the 
LIT method has not been applied to photon scattering, we want to demonstrate with the present work the usefulness 
of this method for this process choosing as a test case the deuteron. That has the advantage that first of all, the 
calculation is quite simple and furthermore allows the comparison with a conventional approach. 

In section II, first we briefly review the formal theory of photon scattering, the low energy limit of the scattering 
amplitude and the concept of generalized polarizabilities as basic quantities. Then we develop the general formalism 
of how these polarizabilities can be calculated by the LIT method. The application of this method on the deuteron 
as a test case is described in section III. For demonstrating the method it suffices to limit the explicit calculation to 
the low energy regime where the El contribution dominates. The results are presented and discussed in section IV 
where we also give a summary and an outlook. 

II. FORMAL DEVELOPMENTS 

We start the formal part with a short resumee of the salient features of photon scattering off a bound many-body 
system (for a more detailed review see e.g. ^^), i.e. we consider the process 

^{k)+N,{P,)^^{k')+Nf{Pf), (1) 

for an incoming photon with momentum k and polarization ex and an outgoing photon with momentum k' and 
polarization e'y while the system makes a transition from an initial state with total momentum Pi and intrinsic state 



\i) to a final state with total momentum Pf with intrinsic state state |/). 



A. The photon scattering amphtude 

In view of the weakness of the electromagnetic interaction one can apply perturbation methods. In the lowest, 
i.e. second order in the e.m. coupling, the scattering amplitude is given by two terms, the contact or two photon 
amplitude (TPA) By\{k',k) and the resonance amplitude (RA) R\'\{k' ,k). A graphical illustration is shown in 
Fig.m 
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FIG. 1: Diagrammatic representation of the resonance (direct (a) and crossed (b)) and the two-photon amplitude (c) for photon 
scattering. 



Accordingly, the total scattering amplitude is the sum of these two contributions 

T[\{k',k)^Bi\ik',k) + R{\ik',k), 
where the two-photon amplitude as depicted in diagram (c) of Fig. [T] has the form 

Bi\{k', k) = -{f\J d^xd?ye^'-^e-^^-ye'^,- B [x, y) ■ ex\i) , 

and the resonance amplitude (RA) (diagrams (a) and (b) of Fig. [1]) is given by 



i?{;,(fc',fc) = (/l 



with the intermediate propagator 



ev • J(-fc', 2Pf + k') G{k + ie) ex ■ J{k, 2P, + k) 

+ex ■ J{k, 2Pf - k) G{-k' + ie) e'^*, ■ J(-fc', 2P, - k') \i) 

G{z) ^ (H - E, - z)-\ 



(2) 



(3) 



(4) 



(5) 



In these expressions, the cm. motion has been separated applying translation invariance. Thus the initial and final 
states refer to the intrinsic motion of the many-particle system only. In principle, the Hamiltonian H and the energy 
Ei contain contributions from the cm. motion. Furthermore, the splitting in a resonance and a two-photon amplitude 
is gauge dependent. This gauge dependence is reflected in the gauge conditions sketched briefly below. 

The cartesian tensor operator B of rank 2 in Eq. ([3| represents the second order term of the e.m. interaction with 
the system under consideration, and the current operator in Eq. ^ 



J(fc,P)=j(fc) + 



2AM 



p{k) 



(6) 



acts on the intrinsic variables of the system only. It consists of the intrinsic current j{k) plus a term taking into 
account the convection current of the separated cm. motion. M denotes the nucleon mass and A the mass number of 
the nucleus. The intrinsic charge and current operators consist of a kinetic or one-body and a meson exchange part 

p{k) = p[i](fc)+p[2](fc), (7) 

J{k) = J[i]{k)+J[2]{k), (8) 

with 

P[i](fc) = 5]e,e'*^S (9) 



'W(^) = ^E (^'iP''^"*'''} + /^"?' X ke-'^-^') . (10) 



2A/ 

I 

Here, e; and fii denote charge and magnetic moment of the l-th particle and pi and ai its internal momentum and spin 
operator. The expressions for the corresponding exchange operators depend on the interaction model. At least in the 
nonrelativistic limit, the exchange contribution to the charge density vanishes (Siegert's hypothesis). Furthermore, 
also the TPA consists of a kinetic one-body contribution and a two-body exchange amplitude 

B{k',k)=B[i]{k',k)+Bi2]ik',k), (11) 

where the kinetic one-body operator is given by 

I 

which is the sum of the individual Thomson amplitudes. 

Gauge invariance of the electromagnetic interaction leads to gauge conditions for the various e.m. operators accord- 
ing to 

k-j{k) = [H,pik)], (13) 

k'-B{k',k) = [p(-fc'),J(fc)], (14) 

where H = T + V denotes the intrinsic Hamiltonian of the nuclear system with T as kinetic energy and the interaction 
potential V. Separating the one-body and exchange contributions, one finds 

fc-J[i](fc) = [r,p[i](fc)], (15) 

fc-J[2](fc) = [F,p[i](fc)] + [r,p[2](fc)], (16) 

k'-Bii]{k',k) = [p[i](-fc'),j[i](fc)], (17) 

k'-B[2]{k',k) = [p[i](-fc'),/[2](^)] + [P[2i(-^'),J[ii(fc)]. (18) 
One important consequence are the low energy limits [ij, [ij] 

J(0) = [H,D], (19) 

S[i],A'A(0,0) - -el'-e^A^, (20) 

S[2],A'a(0,0) = -{i\[e'^1-D,[V,e^-D]]\i), (21) 

NZe'^ 
i?5:,,(0,0) = elT.eA^^-i3f^],,,,(0,0), (22) 

resulting in the low energy limit for the total scattering amplitude 

Tli,(0,0) = -elT-e-A^, (23) 

which is the classical Thomson limit. 



B. Generalized nuclear polarizabilities 

The expansion of the scattering amphtude with respect to the total angular momentum transferred to the nucleus in 
the scattering process leads to the concept of generalized polarizabilities. These polarizabilities allow in a convenient 
manner to separate geometrical aspects related to the angular momentum properties and dynamical effects given by 
the stren gth of the various polarizabilities. To this end, one starts from the multipole expansion of the plane wave 
(see, e.g. [Hi) 

e, e'^- = -V2^Y.^Dij^{R) ^ yAi^iW^- k) , (24) 

LM v=0,l 

with standard electric and magnetic multipole fields A^{M'^; k), where ly indicates the type of multipole field (M*^ = E 
(electric) and M^ = M (magnetic)). Then one expands the current operator in terms of electric (M"'^ = E^) and 
magnetic (M^'-^ = M^) multipole operators 



GA -J 



m = -^/2^Y.^^Mx{R) E ^"^^m"" > (25) 



LM 



where R denotes a rotation which carries the quantization axis into the direction of fc, and D^j^^ denotes the corre- 
sponding rotation matrix [1^. A similar expansion holds for the two-photon operator. 

The electric and magnetic multipole fields of order L of the incoming photon transfer an angular momentum L 
according to the strengths of the corresponding nuclear transition multipole moments. Similarly, the scattered photon 
transfers angular momentum L' . These consecutive momentum transfers can further be classified according to the 
total momentum transfer J to the nucleus with \L — L'\ < J < L + U. The corresponding strength is given by the 
polarizability 

Pl}y\k',k)= J2 y'''yPrf,j{M'''L',M'^L,k',k), (26) 

where v classifies the type of multipole transition as already mentioned. A graphical visualization of the generalized 
polarizability is shown in Fig. [21 




FIG. 2: Graphical representation of the contribution of tlie direct term of the resonance amplitude to the generalized po- 
larizability Pif^j{M'^ L' ,M'^ L,k' ,k) with consecutive angular momentum transfers L and L' by multipoles of type v and v' , 
respectively, coupled to total angular momentum transfer J. 



Then the expansion of the total scattering amplitude in terms of these polarizabilities reads 

// J Ii\ f L L' J 

xPt;f\k\k)Dij,^{R)D^,[,_y{R'), (27) 



r/:\(A?',fc) = (-)i+^'+^/-^^- Y.L',M-,LM,j (-)"+''(2^ + i)(_a'/, i m){i 



where (7^, Mj) and (/y , Afy) refer to the angular momenta and their projections on the quantization axis of the initial 
and final states, respectively. Furthermore, R and R' describe the rotations which carry the quantization axis into the 
directions of the photon momenta k and fc', respectively, and Dj(^j \{R) and Dj(^.p _y{R') the corresponding rotation 
matrices in the convention of Rose [1^] . 

As for the scattering amplitude, the polarizabilities can be separated in a TPA and a resonance contribution 

P,fj{Ar'L', IVrL, k', k) = P^jP/{Ar'L', WL, k', k) + Pl^^jiM"' L' , IvrL, k', k) , (28) 



where for the resonance amphtude one has 

f f' 

P^^'iM" L',M''L,k',k) = 2n{-)^+''^^ 

J 



(29) 



x{IfEf\\ {\^Ar'^^' {k')Gik + ie)Ar'^{k)] + [Ar^^(fc)G(-fc' + ze)iir'^^'(fc')] )\\I^E^) 



Here the symbol "[. . -Y" means that the two multipole operators are coupled to a spherical tensor of rank J. Fur- 
thermore, we have neglected the small cm. current contribution of Eq. ^. The two-photon contribution to the 
polarizability is given by 



/^7/(i^r^^',i^r■^ft^fc) = 2n{-)^+-'+'^{ifEf\\Jd^xd'y 



A^ (iVr ;fc,f)- B (a;,y) • A^(iVr;fc,y) \\hE,) 



(30) 



The evaluation of the TPA contribution to the polarizabilities is straightforward once the TPA operator B {x, y) is 
given. 

For the resonance contribution, one finds by evaluating the reduced matrix element in standard fashion (see e.g. [16|) 

P'7fj(Ar'L',ArL,A;',/c) = 27r(-)^+^^+^'iL' 

L L' J\ {IfEfWM''' '^' {k')\\InEn){InEn\\M''^^{k)\\hE^) 



If li In 



En - E., - k ~ ie 



-(-) 



L+L' + jfL' L J\ {IfEf\\M'^-^{k)\\In){InEn\\Ar'-^'ik')\\l,E,} 



If li In 



E„ - Ei + k' - ie 



(31) 



Obviously, the calculation of the resonance part is more involved because of the summation over all possible interme- 
diate states |/„) and energies En- 

The low energy expansion of the polarizabilities has been discussed in Ref. 14]. For k = Q only the scalar El- 
polarizability is nonvanishing, i.e. 



PjiEl,El) 



fc=0 



^,70 TVs 



e^Z^ 



Ma ' 
with / as ground state spin, which corresponds to the Thomson amplitude. 



(32) 



C. The scattering cross section 



Before turning to the LIT method, we will briefly review the scattering cross section in terms of the polarizabilities. 
For unpolarized photon and target it is given by 



da k' c{k,pi, k') 
'dn ~ T 2(2Ji + 1) 



J\ 2^ KA'A,Af/,Mi(^ '^)l ' 



X,\',Mi,Mf 

with a kinematic factor resulting from the final state phase space and the incoming fiux factor 

c{k,Pi, k') = '-J — . 

In terms of the polarizabilities one finds for the cross section fl3, [la] 
d(T k' c{k,pi, k') 



(33) 



(34) 



dn k 21, 



-^^^^^ J2 Yl P^f,J{M'''L',^'FL)P*fJ{^'FK',N^K)gi^'''^'''^'''^{e), (35) 

'■ T I T Ty I TJ' 1 ,.l . . r-.l r-. 



L'L,K',K,J u' ,i',h'' jU 



where the angular functions are given by 

v'L'vLv'K'vKin\ \)' 



9 J 



m 



(2J + l)(_)i+-K+'''+P' ^(2j + 1)(1 + (_)i+A-+J+^+P)(l + (_)L'+A-'+j+i.'+p'-) 



-T";;)(t^V„){,^-^:5^''^(»-)^ 



with Pj{cos6) as Legendre polynomials. For pure El transitions one obtains 



da{El) k' c{k,pi,k') ^^^_^ ^ e^ 

J (2/. + 1) 1. 1^.^(^1.^1)1 9j (0), 



dn 



where in an abbreviated notation 



(36) 



(37) 



9oiO) = ^(1 + 



6 

gf\e) = ^{2 + sin' 9), 

92\0) = l(13 + cos2. 



(38) 
(39) 
(40) 



D. Application of the Lorentz integral transform method 

A convenient method for the evaluation of the polarizabilities is provided by the Lorentz Integral Transform (LIT) [9| 
as applied to exclusive reactions. For this purpose we separate in Eq. (|29l) the intermediate propagator from the 
reduced matrix element by writing 



/•oo 

G{k + ie)= / dE 
Jeo 



5{H - E) 
E- E,-k-ie 



(41) 



where Eq denotes the ground state energy, and introduce for the separated reduced matrix element as a convenient 
abbreviation a quantity which henceforth will be called the polarizability strength function 



.If I, 



Fl^'L',.L)Ak\KE) = 



_\J+if+h 



(-) 



J 



rv' ,L' /ij 



{IfEfW M" '^ (fc') X S{H - E) KF^^{k) 



\LEi 



(42) 



One should note that in general the strength function is off energy shell, i.e. E ^ Ei + k. Evaluating the reduced 
matrix element, one obtains 



P!^'L',.L)Ak\k,E)^Y.P(^-^FAt /' /} {IfEf\\M'^'-'^'ik')\\I^,E){I^,E\\M'^-^{k)\\I,E.,) 



(43) 



with p{I, E) as density of states for a given energy E and angular momentum /. In terms of the strength functions 
the polarizability becomes 



Pl]'^j{A'r'L',KrL,k',k) ^ 2n{-f+^f+^'LL' 



?ifi^ 



X / dE 

'En 



Ff'''r, r,,(k',k,E) , Ff^!' ,,„j(k,k',E) 



+ {-Y 



E-E,-k-is ' ' E-E, + k' -ie 

One can separate the real and imaginary parts of the propagator according to 

(i^7fj(iir'L',ArL,fc',fc))Re = 2i:{~)^+'f+''LL' 



.Uh 



,Jfii 



xV / dE 
Jeo 



F;V;, r,r(k',k,E) , F^, ,,„j(k,k',E) 



E-E,~k 



E-E, + k' 



(44) 



(45) 



where V stands for the principal value of the integral, and 
{P^f^j{Ar'L',M''L,k',k))i^ = 2^2(_)L+/,+/,^^/ 

X 'P!^'L',.L)Ak', k, E, + k) + (-)^+^'+^4t.'L')./(fc' fc'' ^^ - fc')] • (46) 

where the second term contributes only if Ei > Eq, i.e. if the initial state is an excited state. The subscripts "Re" and 
"Im" indicate the contributions of the real and imaginary parts of the propagator, respectively. For elastic scattering 
the strength function is real and then Eqs. (PS)) and (|46l) represent the real and imaginary parts, respectively, of the 
polarizability. 

The strength functions are the principal quantities which are determined by the Lorentz Integral Transform method. 
Thus the main task is the evaluation of the strength function FJ,^, ^j^\j{k'^ k, E). To this end we first consider the 
following partial strength function for a fixed intermediate total angular momentum state |/„M„) as defined by 

F'j['ilik',k,E) = piI„,E){IfEf\\Ar''^'ik')\\Ir.,E){Ir,,E\\Ar'^ik)\\hE,) . (47) 

In terms of these partial strength functions the polarizability strength is given by 

F;i'L',.L)Ak', k,E)=J2 {/;. t' i} Fl'i^%^^ ^' E) ■ (48) 

The partial strength can be expressed in terms of states of good total angular momentum /„ which are generated by 
the action of a multipole operator M'^^ on a state ip^ with good angular momentum / according to 

liAr'^'ik) X ^')InMn) = lA^r'^^ik) X |/)1 '" . (49) 



Namely, using 

lTl\,T\(T\,JV^L(i,\ ^ ,I^U\j j^,j \ _ ST (_\Ii-L-M„j I 

.M M, -Af, 



(/Af|(Ar^^(fc)xV^')/nAf„) = E(-)''"'^"''"^"(m M- -M ) iIM\Mlt\hM^^ 



6njMM y^\ " (/„||iVr^||/,) , (50) 



one finds 

Fl!i!-:'l(k', k,E) = (^_Y.-i.+L-L'W ^((Ar'^i'(A;) X ^'f)I„AQS{H - E)\{Ar^^{k) x ^^0^«M„) . (51) 

Then we perform a Lorentz integral transform with a complex argument a — <tj^ + iai 

T T T f°° F!jl''''l!Hk',k,E) 

L-'iiiiUM ^ j^^,E^±^^iy-l^ (52) 

Inserting the explicit form for the strength function of Eq. (|5ip and integrating over the (5-function, one finds consec- 
utively 

i'^S(fc',fc,^) - (-)^"-^-+^-^'+^' 

x^((M'^''^'(fc) X ip'i)InhU{H - a)-\H - a*)-\Ar^'^{k) x i>'^)InM^) 

= (-)'""''+''"'''+'^'pan>'^)EW;£M„(fc':'^)i^/;f7„M„(fc:'^)), (53) 

where p{In,(j) takes into account the possibility, that for the given /„ and a several Lorentz states may exist. Here, 
the Lorentz state of good total angular momentum /„ and projection M„ obeys the equation 

(H - a*)\4i-j'Xj,,Jk,a)) = |(M'''^(fc) x ^|/^)In^U} . (54) 

After inversion of the Lorentz transform, one obtains the desired polarizability strength function from Eq. (j48p 
which then serves for the evaluation of the appropriate generalized polarizability. 



III. THE DEUTERON CASE 

For the calculation of deuteron elastic photon scattering the 7-deuteron cm. -system is chosen, where one has k = k' . 
In order to demonstrate the LIT method it suffices to consider the at low energies dominant El transitions with the 
El operator in the long wave length approximation (Siegert form) 

El, - i[H, D\,] , where D\, = 0= r Y^t {^) (55) 

is independent of k. Here a denotes the fine structure constant, and (r, fi) represents the relative neutron-proton 
coordinate. Furthermore, H denotes the intrinsic two-nucleon hamiltonian containing the nucleon-nucleon interaction. 
The small cm. contribution to the Hamiltonian is neglected for simplicity. Due to the dipole approximation only the 
polarizabilities Pj{El,k) (J = 0,1,2, in an abbreviated notation) contribute. Instead of the corresponding partial 
i?l-strength function 

Fh^^EiiE) = i-y-'J2((^D\H]xi;',)jm\6{H-E)\i[D\H]xi:',)jm) 

= {-y-\E Eof J2((^' ^ ^^)^'"^l ^(^ ^) l(^' ^ ^^)^'"^) ' (56) 

ni 

where Eq denotes the ground state energy, we will consider the reduced partial strength function 

One should note that F^-^^-^^-^ is independent of k. The associated Lorentz state obeys as LIT equation (see Eq. (15^ ) 

{H - <J*)\i^jm{a)) =\{D^x ^i)jm) , (58) 

which is independent of k. The Lorentz state can be expanded into partial waves according to the orbital angular 
momentum I 

(r,f]|V',™(a)) = ^ E $^,(a,r)(0|(/l)jm), (59) 

where | (n)j to) represents the spin-angular state of orbital angular momentum I coupled with total spin-one to a total 
angular momentum j, and the matrix element refers to spin and angular degrees of freedom only. The state ^jm(c) 
generates the LIT L^{-'^^ of the reduced strength F^{-'^-^. 

Inserting the expansion of ^jm into Eq. (|58p and projecting onto a state | (/l)j to,) , one finds a set of radial differential 
equations 

I (^ - ^) - -*] ^A-.r)+'ZV,,,,^,,{-,r) = ^^rf,,{r) (60) 

with 



f,i{r)^6nu{r)+ (-)^+i 3 v^f (^^ 1 J^ |2 1 jj ^^^^ ^ ^g^^ 

where uir) and w{r) are the deuteron radial s- and d-wave functions. Two of them are uncoupled (^Pq '■ ^ — ^tJ — ^') 
•^Pi : I = l,j = 1) and one is coupled (^^^2 — ^-f2 : h = l^h = 3,j ~ 2). Equation (IMl) is very similar to a radial 
Schrodinger equation, but with a complex energy a and an additional source term on the right-hand side. 
For convenience we introduce the following three reduced LITs 



L,{a) := (-)^-i 4i-f^^(a) ^ ^(^^„(a)|^,„,(a)) = a^ / \^ ,,{a , r)\'' dr ^ 



(62) 



and corresponding reduced strength Fj{E) — F^^-'^j^{E)/(2j + 1). Then the polarizabiUty strength function becomes 



^"- ^--^^^^eh^'^^(2j-+i)II \ i\p^(E) 



FEl^EliE) 



An 



1 1 J 



(63) 



and, according to Eq. p6)) . the El polarizabihties are given by 



{p:r{Ei,k))i 



= -67:' F^'-^^^.ik + Eo) 



iPr{El,k))ne = ^^(-)^(2j + l)^J I '^^\V jdE{E-E,YF,{E) 



E — En 



(64) 



{-Y 



E-En 



= -7^ *'(Pr(i^M')W,,,_,.,,^, 



(65) 



The latter expression for the real part corresponds to the dispersion theoretic approach of Ref. [1| . It is a consequence 
of the fact, that in this special case of taking the _El-operator in the low energy limit the polarizability strength 
EeYei becomes independent of the photon momentum k. However, in general this is not true. 



IV. 



RESULTS AND DISCUSSION 



For the numerical solution of the radial equation (pO]) for the radial Lorentz states ^ji {a, r) we have chosen the 
Argonne potential AV18 [l^ as interaction model. The resulting three Lj{a) are shown in Fig. [3] for a constant 
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FIG. 3: (Color online) Lorentz integral transforms Lj{a) with a/ = 5 MeV: j — (dotted), j = 1 (solid), J = 2 (dashed). 



CT/ = 5 MeV as function of aji up to 100 MeV. The figure shows that the three LITs have quite a similar behavior. All 
three exhibit a pronounced peak in the low a^ region, only the peak heights are slightly different. For the principal 
value integral in (I45p also high-energy contributions could play an important role, therefore we illustrate in Fig. 2] the 
transforms in an extended an range. One finds that the behavior at high a^j is approximately described by a~^ . It 
shows that at large an the transforms are dominated by low-energy contributions. In fact, for the extreme case of a 
(5-strength, i.e. F{E) = Fq S{E - Eq), one obtains L{a) = Fo/{{aR - EqY + crj). 

One also notices in Fig. 2] the onset of oscillations at higher an, which are more pronounced for Li and L2- The 
origin of these oscillations lies in the relatively small value of 5 MeV for aj, which makes a high-precision solution of 
(|60p with increasing a^t more and more difficult. Such a small ai value, however, is advantegeous in the low a^ region, 
where (i) it does not lead to numerical problems and (ii) prominent structures of a small width could be present, e.g. 
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FIG. 4: (Color online) As in Fig. |3l but for an extended range of an. 



a low-energy peak. The small ct/ value will then allow one to resolve such details and nonetheless will not lead to 
problems for the reconstruction of the high-energy strength, since structures with a small width are not present at 
higher energies. One could also completely avoid the oscillations at higher an by choosing a transform with different 
CT/ values for low- and high-energy regions, as has been done in [20[. 

In order to obtain the strength function F^jE) one has to invert the integral transforms Lj defined in ((62|) . Details 



about the inversion of the LIT are found in [9|, |2l|, and further, more general inversion aspects are discussed in [2 

Accordingly, wc use expansions of the calculated LlTs, Lj{aR,(7j ~ 5McV), in a set of basis functions x„ {n = 
1, 2, ..., N), where the expansion coefficients arc determined by a best fit. Here we take, as in Q, the set (for fixed ct/) 



xi'K^R) 



dE- 



x'n\E) 



{E - anf + aj 



with 



x'^^\E)^E'-'exp{^ 



a2E^ 



(66) 



(67) 



where a^ and /3 are nonlinear parameters. The inversion results, Fj(E), are shown in Fig.[5l One observes pronounced 
low-energy peaks with a strong subsequent fall-off which becomes weaker at somewhat higher energies. The fall-off 
becomes stronger again at even higher energies, namely beyond about 200 MeV (300 MeV) in case of Fq (F1/2). As a 
matter of fact, the inversion results are already not very good somewhat below those energies, since there we do not 
find, as it would be necessary, a stability for the inversion solutions for a limited range of the number of basis functions 
TV. The origin of this problematic high-energy behavior lies in the choice of basis functions. By construction they all 
have an exponential high-energy fall-off and thus are not suitable to describe a function with a different high-energy 
behavior, since N cannot be increased arbitrarily because of the numerical accuracy of the transform to be fitted. 

As already mentioned, for the principal value integral in (|45p also high-energy contributions of the strength function 
could matter. Therefore it is better to search for a basis set which does not have the shortcomings of the set in (|67p . 
but which is more appropriate to describe the strength function over a much larger energy range, even if it could lead 
to a somewhat less precise inversion at lower energies. For this purpose we introduce an alternative set without an 
exponential fall-off: 



An -^ 



with 



{E) = {n + -)+[(3-in+-) 



E 



(68) 



(69) 



where /3 < 0, 7, and i?asy are nonlinear parameters. In the present case we have calculated the LIT up to 1000 MeV 
and thus we set £^asy = 1000 MeV, i.e. all basis functions have an asymptotic fall-off with E^ at E = 1000 MeV. The 
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FIG. 5: (Color online) Strength function Fj{E) for j = (dotted), j = 1 (solid), j = 2 (dashed). 
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FIG. 6: (Color online) Strength functions Fj from inversions with first (thin lines) and second basis set (thick lines): j 
(dotted), j = 1 (solid), j = 2 (dashed). 



inversions with basis set {x*-^''} are shown in Fig. [6l At lower energies one finds a similar picture as in Fig. [5l the 
high-energy behavior, however, is quite different. One further notes that a constant asymptotic high-energy fall-off is 
already reached at energies considerably lower than iJasy = 1000 MeV showing that this set is a very good choice. 

Figure [6] also shows that at lower energies the inversion results do not depend significantly on the choice of the 
basis set. On the other hand the precision of the agreement is difficult to judge on a logarithmic scale. Therefore, we 
show in Fig. [7] the various inversions at lower energies in a more detailed formn, i.e. as relative differences 



A,iE) 



(fI^'He) 



FPiE))/FP{E). 



(70) 



where F- and f!- correspond to the inversions with the first and the second basis set, respectively. At very low 
energies one finds the largest differences for Fi, namely about 7%. This relatively large difference is due to the steep 
rise of the strength function right above threshold. Here we would like to mention that, as suggested before, our first 
basis set leads to a better fit of the calculated LITs in the low ajj, region. Furthermore, Fig. [7] shows that up to 50 
MeV further differences remain very small, in fact less than about 1%. For Fi and F2 the picture does not change 
very much for higher energies, whereas in case of Fq the difference increases quite substantially beyond 80 MeV. This 
can be interpreted as a precursor effect for the strong fall-off beyond 200 MeV which leads to a somewhat oscillating 
inversion result already at considerably lower energies. In order to have the best description in the whole energy range 
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FIG. 7: (Color online) Relative differences Aj{E) of the strength functions shown in Fig. |6l Notation of curves as in Fig. [G] 



we combine the inversion results of the two basis sets by taking F (E) for E < Ei and F- (E) for E > Ei with Ei 
equal to 30, 18, and 33 MeV for j=0, 1, 2, respectively. 

The strength functions Fj can also be obtained from a standard calculation of deuteron photodisintegration. Ac- 
cording to Ref. [1[ one has 



linPj{El,El) = 47rfcV(-) 



j 



1 1 J 
1 1 J 



E 

fils 



E\fi,jls)\' 



(71) 



where the Sl-matrix elements are defined by the multipolc expansion of the total unpolarized photo absorption cross 
section [H, [11 



CTo = 



(47r)2 



E 



1 



3 ^^ 2i + l 



{\E^{f^,jls)\' + \M^{f^,jls)\'). 



(72) 



Here, the quantum numbers {fJ.ls) classify the various components of a two-body scattering solution with total angular 
momentum j j23l . [2J] . Comparison with Eq. (|64p yields the relation 



Fj{k + Eo) 



8 



3fc(2i + 1; 



j2\E\f^,jism\' 

flls 



(73) 



A comparison of the standard calculation with the LIT approach is shown in Fig.[8l The agreement is quite satisfactory 
in view of the fact, that in the standard calculation the complete £^l-operator is used and not its low energy form as 
in the present LIT approach. 

Having obtained the reduced strength Fj{E), one can then determine with the help of Eqs. (|M)) and (|55|) the 
imaginary and real parts of the polarizabilities PJ°^(_E1, El, k). In Fig.lHlwe compare our results for the polarizabilities 
to those of [IJ. One should note that the real parts are normalized to zero at fc = for J = and 2, which is not 
necessary for J = 1 since one has Pj^^{El,El,k = 0) = 0. This normalization takes into account implicitly the 
neglected contributions of the two-photon amplitude at fc = 0, which are needed in order to comply with the low 
energy theorem of Eq. (j32]) for J — 2, whereas for the scalar polarizability we have to add the classical Thomson limit 
to the normalized resonance part according to 



3e^ 

Re (PoiEl, El, k)) = Re {Pn^^El, El, fc)) 



(74) 



where md is the deuteron mass. This procedure is justified according to the discussion in [1^ (see Fig. 3 of [18j). 
where it is shown that further contributions to the two-photon amplitude beyond the low energy limit are negligibly 
small for fc < 40 MeV and remain quite small up to about 60 MeV. Thus, for calculations below 60 MeV, it seems to 
be quite safe to simply add the Thomson term to the normalized resonance ElEl scalar polarizability. 
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FIG. 8: (Color online) Strengths functions Fj{E) weighted with {E — Eo) for the Argonne V18 potential in standard approach 
with complete _El-operator (thin lines) and LIT result with low energy _El-operator (thick lines): j = (dotted), j — 1 (solid), 
j = 2 (dashed). 



With respect to the comparison in Fig. |9l we should mention that the polarizabilities of [1| are calculated in the 
Breit system. The relation between the photon momentum k of the Breit to that of the cm system depends on the 
scattering angle. Only in forward direction they turn out to be exactly the same. In addition it should be mentioned 
that in [1] a more complete calculation has been made where the full nuclear one-body current and also pion exchange 
currents have been taken into account. Furthermore, in [ij a different NN potential has been used. Despite these 
differences, one finds a rather good agreement between the two calculations for the scalar and tensorpolarizabilities. 
Larger differences are found for the vector polarizability, for which real and imaginary parts from jl| are somewhat 
larger than our results. In Fig. 9 also the imaginary parts of the Pj evaluated by the standard method are shown. 
For J — and 2 they are very similar to our results, actually the agreement is even better than one would expect 
from the results of Fig. 8. On the other hand, for J = 1 the standard result agrees quite well with that of Ref. [l|. 
However, one should keep in mind, that the vector polarizability is about an order of magnitude smaller than the 
scalar one due to large cancellations between the Fj . Thus the small differences in the Fj according to Fig. |5] result 
in considerably larger differences in Pi. 

With the calculated polarizabilities we can determine the cross section for pure El transitions given in (j37p . In 
Fig. [TU] we the show the resulting cross section in comparison to experimental data at fciab = 55 MeV. One sees that 
our theoretical result overestimates slightly the experimental cross section, although within the experimental errors. 
Since data exist only for a very limited angular range one cannot draw any conclusion whether the ElEl cross section 
reproduces the correct angular shape. However, from other theoretical calculations 1, 3-6] mentioned above, it is 
known that additional polarizabilities PY^{M^ L',M'^L, k) from neglected other multipoles, e.g. Ml and i?2, play a 
considerably less important role at lower energies. Furthermore, in Refs. 0-3 , it is pointed out that also the internal 
nucleon polarizabilities lead to non-negligible contributions. In fact, one aim of present-day Compton scattering 
experiments on light nuclei is the determination of the nucleon static electric and magnetic dipole polarizabilities, in 
particular the ones of the neutron. 

We summarize our work as follows. After a brief overview over the theory of photon scattering, we have described as 
the central issue of our work the application of the LIT method to the calculation of the photon scattering amplitude. 
This method has the great advantage over conventional methods, since the resonance amplitude can be evaluated in a 
very efficient way, in which wave functions of the continuum spectrum of a particle system have not to be calculated 
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FIG. 9: (Color online) Imaginary (left panels) and real (right panels) parts of the polarizabilities P}''^{El,El,k) (solid) in 
comparison to those of 1] (dashed): J — (top panels), J — 1 (middle panels), J — 2 (bottom panels). Note that the real 
parts are normalized to zero at fc = (see text). For the imaginary parts also the results with from the standard calculation 
are shown (dotted) 
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FIG. 10: Comparison of the differential scattering cross section in the unretarded dipole approximation with experimental data 
from I^ (fcub = 55 MeV). 



explicitly. Because of this reason, and in contrast to conventional state-of-the-art methods, the LIT approach can be 
applied also to systems with more than three particles. 

As a first application and test case, we have considered elastic photon scattering off the deuteron at low energies. 
For this case we have taken into account only El transitions in the unretarded dipole approximation. In order to apply 
the LIT method one has to invert the calculated integral transforms of the so-called strength functions. This inversion 
is discussed in great detail. In particular the high-energy behavior of the strength functions is studied with great 
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care. We have compared the resulting three ElEl polarizabihtics with results from conventional calculations. The 
comparison has shown the reliability of the obtained results. At last we have confronted the ElEl cross section with 
experimental data at a photon energy of 55 MeV. In view of the approximations of the present work the agreement 
is quite satisfactory. This is certainly a strong encouragement to consider in future work photon scattering off more 
complex light nuclei. 
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